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Abstract 



Oh 

O ! Analytical analysis of spatially extended autocatalytic and hypercyclic systems is 

^ ' presented. It is shown that spatially explicit systems in the form of reaction-diffusion 

equations with global regulation possess the same major qualitative features as 
the corresponding local models. In particular, using the introduced notion of the 
stability in the mean integral sense we prove the competitive exclusion principle for 
^ ■ the autocatalytic system and the permanence for the hypercycle system. Existence 

, and stability of stationary solutions are studied. For some parameter values it is 

\ proved that stable spatially non-uniform solutions appear. 

rn , 

Keywords: Autocatalytic system, hypercycle, reaction-diffusion, non-uniform 
^ . stationary solutions, stability 

o 

> 

X 

■ In 1971 Manfred Eigen published a seminal paper on the evolution of error-prone self- 
replicating macromolecules His theory was expanded significantly later on, primarily 
in works of Eigen, Schuster and co-workers [ini IHl US]. One of the principal findings 
was the existence of the error threshold^ i.e., the critical mutation rate such that the 
equilibrium population of macromolecules (the quasispecies in the terminology of Eigen et 
al.) cannot provide conditions for evolution if the fidelity of copying falls below this critical 
level. This critical mutation rate depends on the length of macromolecules and hence puts 
limits on the amount of information that can be carried by a given macromolecule. To 
improve fidelity one needs longer sequences (e.g., a more efficient replicase), to have longer 
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Figure 1: The hypercycle [12]. Each macromolecule (Mj) helps to replicate another one, 
Mj+i, M„ macromolecule promotes the replication of Mi closing the loop; ai, i = 1, . . . ,n 
denote reaction rates. 

sequences one needs better fidelity, hence the chicken-egg problem. An easy and obvious 
solution to this problem is that the early primordial genomes must have consisted of 
independently replicating entities, which, generally speaking, would compete with each 
other (see, e.g., |2Tj and references therein). 

If we consider a simple mathematical description of independent competing replicators 
then the usual differential equations for the growth take the following form: 

— = Qivf - fi{t), i = l,...,n, (1.1) 

where Vi = Vi{t) is the concentration of the i-th type of macromolecules, is the rate of 
replicating, p > is the degree of auto-catalysis, and /i(t) is the term which is necessary 
to keep the total concentration constant, this term depends only on t and not on the 
index, /i (t) = Y17=i (^i'^i'^i the present case; easy to see that this is equivalent to the 
condition Yl^=i "^i = 1- the case if p 7^ 1 we have system with non-linear growth rates, 
which model different coupling strength of the various components, for the discussion of 
such growth rates see, e.g., [231 [21|- Hereinbelow we consider mainly p > (or even, 
p > 1) but remark that p = gives the exponential growth, p = 1 gives the standard 
hyperbolic growth (autocatalysis), and for p < the parabolic growth occurs [27\. It 
is straightforward to show that for p > only one replicator present at t — >■ cxo, the 
competition winds up in the competitive exclusion of all but one types, i.e., the genome 
composed of independently replicating entities is not vital. 

To resolve this situation Eigen and Schuster [12] suggested a concept of the hypercycle, 
a group of self-replicating macromolecules that catalyze each other in a cyclic manner: 
the first type helps the second one, the second type helps the third, etc, and the last 
type helps the first one closing the loop (see Fig. [T]). An analogue to system fll.ip can be 
written in the form 

— = Oi<_i - /2(i), i = l,...,n, (1.2) 
where index coincides with n, f2{t) = J27=i (^i'^f-i'^i- For p = 1 we obtain the standard 
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hypercycle model [IH]. It is known that fll.2l) is permanent, i.e., all the concentrations are 
separated from zero, and hence different replicators coexist in this model. More exactly, 
for short hypercycles, n = 2,3,4, the internal equilibrium is globally stable, for longer 
hypercycles, n > 4, a globally stable limit cycle appears [T7] . 

The problem with the hypercycle model (11. 2p is its vulnerability to the invasion of 
parasites [22] ■ 

We remark that models f 1 1.1 1) and (11.21) are systems of ordinary differential equation 
(ODEs), i.e., they are mean-field models. As a solution to the parasite invasion problem 
it was suggested that heterogeneous population structure can strengthen persistence of 
the system. One of the suggested solution was spatially explicit models [UEllTj, see also 
[21 E] for reviews of the pertinent work. Two major approaches to spatially explicit mod- 
els are react ion- diffusion equations and cellular automata models, and they both were 
considered in the cited works. Which was lacking, however, is an analytical treatment 
of the resulting systems, because in both cases the researchers have resorted to extensive 
numerical simulations. An only notable exception to our knowledge is [32] , where some of 
the models with explicit space are analyzed analytically. An interest in cluster-like solu- 
tions of reaction-diffusion systems resulted in the analysis of spatially explicit hypercycle 
in infinite space [301 EI] • 

Note that models (II. ip and (ll.2p are a special case of the general replicator equation 
|19] . for which several approaches are known to incorporate an explicit spatial structure, 
albeit there is no universally accepted way of incorporating dispersal effects. The solution 
to the problem with equal diffusion rates is straightforward, in this case we, following 
ecological approach, can just add the Laplace operator to the right hand sides of (II. ID 
or (II. 2p . This was used, e.g., in the classical paper by Fisher pi] to model the effect of 
the spatial structure on the invasion properties of an advantageous gene; this approach 
later was generalized by Hadeler p^. However, for the primordial world, it would be a 
too stringent an assumption to have all the diffusion coefficients equal. To overcome this 
difficulty, Vickers et al. introduced a special form of the population regulation to allow for 
different diffusion rates [5l[20l[29j, now in the subject area of evolutionary game dynamics. 
In these works a nonlinear term is used that provides local regulation of the populations 
under question, although no particular biological mechanism is known that lets individuals 
adapt their per capita birth and death rates to local circumstances [13]. In our view, it 
is more natural to assume the global regulation of the populations, hence following along 
the lines of thought that brought to the models (II. ID and (II. 2p . Mathematically it means 
that we assume that the total populations satisfy the following condition 



where x G is a spatial variable now. This approach was first used in |32j. Which 

is important here is that this approach allows to obtain some analytical insights of the 
systems [1]. 

In this text our goal is to present an analytical treatment of the models of prebiotic 
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macromolecules with self- and hypercyclic catalysis with an explicit spatial structure and 
global population regulation in the form of reaction-diffusion equations. 



2 The mathematical models 

Let fl he a. bounded domain, Q C M™, m = 1,2,3, with a piecewise-smooth boundary F. 
The spatially explicit analogue to (11. II) is given by the following reaction-diffusion system 

dtVi = Vi{aivf - fi(t)) + diAvi, i = l,...,n, t > 0. (2.1) 

Here Vi = Vi{x,t), x G fi, t > 0, 9^ = ^, A is the Laplace operator, in the Cartesian 
coordinates A = The initial conditions are Vi{x,0) = V5j(x), p > (although 

we note, that in each particular case we shall specify admissible values of p), and the form 
of fi{t) will be determined later. 

A slight modification of (12. ip gives the hypercyclic system 

dtVi = Vi{aivf_^ - /2(t)) + diAvi, i = l,...,n, t > 0, (2.2) 

where Vq = f„. 

In both problems (12. ip and (12.21) the functions Vi{x, t) are assumed to be nonnegative, 
since they represent relative concentrations of different macromolecules. 

It is natural to assume that we consider closed systems (see also [32]), i.e., we have 
the boundary conditions 

p- =0, z = l,...,n, (2.3) 

where n is the normal vector to the boundary F. 

It is assumed that the global regulation of the total concentration of macromolecules 
occurs in the system such that 



n „ 

/ Vi{x, t) dx 

i=l 



(2.4) 



for any time moment t. This condition is an analogous condition for the total concentra- 
tion of replicators in the finite-dimensional case [TH] . From the boundary condition (12.31) 
and the integral invariant (12.40 the expressions for the functions /i(t) and /2(i) follow: 



n „ 



vr {x,t)dx (2.5) 
and 

n „ 

/2(^) = / aiv'^^_^{x,t)vi{x,t)dx. (2.6) 
i=i 
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Finally we have a mixed problem for a system of semilinear parabolic equations with the 
integral invariant (12 ■4p and functionals fl2.5l) and fl2.6p . 

Suppose that for any fixed x ^ Q each function Vi{x,t) is differentiable with respect 
to variable t, and belongs to the space Hp_^_i{Q) as the function of x for any fixed t > 0. 
Here HL^ is the space of functions with the norm 



\u{x)\f+Ux 



p+1 



+ 



E 

k=l 



du 



dxk 



dx 



Note that if p > 1 then Hp_^-^^{Q) C HK^I), where H2{^1) is the Sobolev space of square- 
integrable functions for which their first partial derivatives are also square- integrable [25] . 

Without loss of generality we shall assume further that volume of the domain Q is 
equal 1: |f2| = 1. 

Our main goal is to analyze existence and stability of the steady state solutions to 
f l2.ip and 02.20 . The steady state solutions are given by the solutions to the following 
elliptic problems: 

diAui + Ui{aiU^^ - fi) = 0, i = l,...,n, (2.7) 

and 

djAuj + Mi(ajMf_i - /s) = 0, i = l,...,n, mq = u„, (2.8) 

with the boundary conditions d^Ui = on F; Ui{x) G Hp^^{Q). The integral invariant 
(12. 4p now reads 



/ Ui{x) dx = 1, 

i=l 



(2.9) 



the values of /i and /2 are constant: 

n 

and 



i=l 



aiU\ 



p+1, 



x) dx 



f2 = y2 a'iU^-iix)ui{x) dx. 
i=i 



(2.10) 



(2.11) 



If it is assumed that di = d2 = ■ ■ ■ = d^ = ^ then the equilibrium points of (11.11) 
and (II. 2p coincide with the steady state solutions to (12.10 and (12. 2p . These solutions are 
spatially homogeneous. The converse is also true: the spatially homogeneous equilibria 
of systems (12.10 and (12. 2p are fixed points of the dynamical systems (II. ip and (11.20 
respectively. 

The coordinates of these spatially homogeneous solutions are straightforward to write 
down. Let Pi = {ai)~p and consider the sum /? = ^/Sj, where the index of summation is 
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determined later. All spatially homogeneous solutions to (12.71) are given by 

1 

Qjk = ^(/^i' • • • ' 0' /^j+i! Pk-i, 0, /5fc+i, . . . , 



ending with the vertices i?j = (0, . . . , 0, 1, 0, . . . , 0) (unity at the z-th place) of the simplex 
Yl^=i = 1) ^ 0; for each steady state (3 in obtained by summing through all non-zero 
elements in the vector. 

The spatially homogeneous stationary solution to (12.81) is given by 

3 Stability of spatially homogeneous equilibria 

Let M° = {ui, . . . be a spatially homogeneous solution to system (12. ip . In the usual 
way we assume that the Cauchy data are perturbed 

(Pi{x) = u° + Wi{x), i = l,...,n. 

Here Wi{x) G Hp_^_^{fl). Inasmuch as we have 

1=1 

then from (12. 4p it follows that 

n „ 

y / Wi{x) dx = 0. (3.1) 
i=i 

Consider the following eigenvalue problem 

A^(x) + \^{x) = 0, X G 9n^(x)|:,er = 0. (3.2) 

The system of eigenfunctions of this problem ipQ^x) = 1, {ipi^x)}'^^ forms a complete 
system in the Sobolev space Hl{Q) [25] such that 

{i!i{x),iljj{x)) = / ipi{x)iljj{x)dx = 6ij, 
'n 
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where 6ij is the Kronecker symbol. The corresponding eigenvalues satisfy the condition 
= Ao < Ai < A2 < . . . < Aj < . . . , lim Aj = 00. 

Hence for p > 1 we assume that Wi{x) can be represented as 

00 

= '^c](l)j{x), i = l,...,n, (3.3) 

j=0 

where c* are constant. 

Denote Hs the set of functions w{x) G Hl{Q) such that < S, where 6 > 0. 

Theorem 3.1. For p > 1 all spatially homogeneous stationary solutions to (12.11) are 
unstable with respect to any perturbation from the set Hs if 

0< — <f-, i = l,...,n. (3.4) 

Qi Ai 

The solutions Ri = (0, . . . , 0, 1, 0, . . . , 0), i = 1, . . . ,n are stable when 

^>f. (3.5) 

tti Ai 

Here Ai is the first non-zero eigenvalue of the problem (13.21) . 

Proof. Let W{x,t) = (Wi{x,t), . . . ,Wn{x,t)) be a vector-function belonging to Hg for 
any fixed t. Using (13.11) and (13. 2p we can seek the solution to (12. ip in the following form: 

n 

Vi{x,t) = u'i+Wi{x,t), Wi{x,t) = Y,c]{t)^j{x). (3.6) 

j=0 

Substituting (13.60 into (12.10 and retaining in the usual way only linear terms with respect 
to Wi we obtain the following equations: 

dtWi = ip+ l)ai{u^YW^ - hWi + di/\Wi, W,{x, 0) = ipi{x) G Hs, d^Wi = on T. 

(3.7) 

Consider first the case = Pi. Direct calculations show that fi = {3 ^. 
Multiplying equations (13.70 one after another by the functions ipj and integrating with 
respect to x G f2 we obtain the following system of ordinary differential equations: 

-^ = c]{t){a,pl3-P-dA,), z = l,...,n, j = 0,1,2,... (3.8) 
For j = one has 

= 4(0)exp(a,p/3-Pt), 
therefore, Co(t) ^ cxd as t — 00, which implies that Pi is unstable. 
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Using the same approach it is straightforward to show that Qj, Qjk, ■ ■ ■ are also un- 
stable. 

Now we deal with First note that from (12.91) it follows that 



J2cm = 0. (3.9) 



d{t){ai - diXj), forj = z. 



For Ri we have 

dciit) 
dt 

and 

dciit) 
dt 

Therefore, for j 7^ i, c](t) ^ when t —>■ 00. Taking into account (13.91) we obtain that 
4{t) ^0. If i = j and (IS3D holds then c}(t) 00, if holds then c}(t) 0, which 
proves the theorem. □ 

Theorem 3.2. Ifp > 1 then spatially homogeneous stationary solution P2 to system (12.21) 
is unstable with respect to any perturbations from the set Hs when 

Proof. As before we will look for a solution to (12. 2p in the form (13.61) . After substituting 
(13. 6p into (12. 2p . multiplying by ipj and integrating, we obtain the following system of 
ordinary differential equations for c*(t): 

dc^it) D a,- 

= ^^dr\t) - dA,dj{t), t = l,...n,n + l = l,0 = n, j = 0, 1, 2, . . . 

(3.11) 

Applying the Routh-Hurwitz criterion we obtain that the solutions to (13. lip go to 00 if 
(I3.10p holds, which implies instability of P2- CH 

Remark 3.1. Inverse inequality to (13.101) provides stability of P2 only in the cases 
n = 2, 3, 4. Actually, for j = we have that (13. lip takes the form 

dc^tl _ ^^ r^U . . _ . 

All eigenvalues can be easily evaluated because the corresponding matrix is circular: 

l^j = ^Pj^ j = 0,...,n-l, 

where pj is the j-th root of the equation = 1. The eigenvector (1, 1, ... , 1) does not 
satisfy (13.91) . therefore we exclude it from the consideration. When n = 2, 3 all eigenvalues 
have negative real part, in the case n = 4 P2 also will be stable jT^. For n > 5 there is 
at least one eigenvalue with positive real part, which proves the claim that P2 is unstable 
when n > 5. 
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4 Existence of spatially nonuniform stationary solu- 
tions to systems ( 12. 11) and ( 12.21 ) in one-dimensional 
case 

Here we will prove that when the space is one dimensional, = [0, 1], A = dx-, the models 
(12. ip and f l2.2p possess non-uniform stationary solutions under some additional conditions. 
The boundary conditions now take the form dxVi{Q,t) = dxVi{l,t) = 0. 

Theorem 4.1. For < p < 2 a spatially non-uniform stationary solution to (12. ip exists 
if the following inequality holds 



Proof. We start the proof noting that the dependence of the concentrations Ui{x) in 
f l2.7p on other concentrations and their total regulations occur only through the integral 
invariant fl2.1Up . which does not depend on x. Therefore we can assume without loss of 
generality that each Ui depends on its own variable Xi G [0, 1]. Hence we rewrite ( 12.90 and 
(l2:Tn|l in the form 



/ Ui{Xi)dx, 

i=l 

n „ 

fi = y2 aiU^i^^{xi)dxi. 
i=i 



1, 

(4.2) 



Each equation of system (12.70 can be put in the following form: 

d^r dS^' 

System (14. 3p is a Hamiltonian system for any i = 1, . . . ,n, in which Xi is considered as a 
"time" variable, with the Hamiltonian 



2 di \p + 2 ' 2 
The phase orbits of (14.30 can be found from the standard formula 



V, = ±J2{Hf-Ui), 
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Figure 2: The phase portrait of system (14.31) . The closed curves surrounding A correspond 
to the candidates for spatially non-homogeneous stationary solutions to (12.11) 



where 




From the form of the phase orbits (see Fig. [2]) it immediately follows that there exist 
orbits that satisfy the condition 



Fig. [21 Different diffusion coefficients correspond to the motion along the phase orbits 
with different velocities. 

To prove the theorem we need to show that there exist two values x] and xf such 
that \x} — x^|=l, for i = 2, . . . ,n, and corresponding solutions to (12. 7p satisfy the first 
condition in (14. 2p . 

The solutions to system (12.71) can be found in the explicit parametric form |26j : 






Ui{Xi) 



p + 2 
2ai 



fi r, r>0 



P 




■T 



dt 



(4.4) 



p + 2 
2ai 



fi r, r < 0, 



Ui{Xi) 
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To proceed we need the following lemma (the proof is given in the Appendix). 

Lemma 4.1. The equation Pi{t) = cj+t^—t^^'^ = has two real positive roots < t} < rf 
for all values 

2 



2 + p 







Moreover, 



2 + p 



2 + p 



1 



(4.5) 



(4.6) 



An analogous lemma holds for Qi{t). 

Now we return to the parametric representation (14. 4p . Consider the first derivative of 
the functions m,-: 



dui 
dxi 



dui dr 
dr dxj 



'p + 2- 
2ai 



'd 



Jl 



This expression vanishes at the points tI and rf. Using (14.41) for Xi = and Xj = 1, we 
obtain 



Xi 



Xi 



1 



di 

fi 

di 
fi 



dt 



dt 



+ c 



(4.7) 



1. 



Letting tq = r/ we obtain that c] 

1} 



VW) 

0. We will use the following notation: 
dt fK 



(4i 



Remark that these integrals will exist because the roots of Pi{t) are simple when cj satisfy 
(14. 5p . The formula (14.81) establishes the connection between the values of the constant cj 
and values of the diffusion coefficient di, which determines the velocity of motion of phase 
points. The latter implies that (14. 8 p guaranteers that the motion from the initial point 
{uj, 0) to the final point {uf, 0) occurs during the unit time. 

Now we are going to prove that the solution (14. 4p satisfies the first condition in (14. 2p : 



y2 / Ui{xi)dx. 
i=i -^0 

From ( 149|) it follows that 
'p + 2 



P 



i=l 



h 

ai 



tdt 



1. 



(4.9) 



i=l 



ai 



'di {p + 2)tP+^ /2 - tdt 
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Indeed, we have 



{p + 2)tP+V2 - tdt _ 1 f^' d{tP+^ - 
I VW) "2 

1 r" d{p,{t)) 



Pdrf 



Pr(r} 



0. 



2 7.1 

Using (14.81) to find /i and substituting this expression into (14.91) we obtain 

> + 2" 



E 

i=l 



(4.10) 



where 



dt 



tdt 



To conclude the proof we need the following lemma, the proof of which is given in the 
Appendix 

Lemma 4.2. If0<p<2 then the following inequality holds: 



p + 2 



I > 



TT 



P 



{4.11] 



Applying the result of Lemma 14.21 to (I4.10p we obtain that if (14. ip holds then there 
exists a spatially non- uniform stationary solution to (12. ip . □ 



Remark 4.1. From the symmetry of the system, the time needed to get from the point 
{u},0) to the point (m^,0) is the same as the time needed to get from (^^,0) to (^^,0), 
and the speed of movement is inversely proportional to y/d^. Therefore, reducing these 
values twice we guaranteer that spatially non-uniform stationary solution exists, which 
corresponds to the full cycle in the phase plane; reducing 4 times we obtain the solution 
which corresponds to the movement of the phase point along the cycle two times, and so 
on. Hence system (12.11) has non-uniform stationary solutions that correspond to movement 
along the cycles in Fig. [2] arbitrary number of times (see Fig. 

Remark 4.2. We introduce the following parameter 

1=1 ^ ' 

Theorem 14.11 can be restated as follows: If 

d < jji = (p7r~^)p 



12 



0.2 0.4 Xi 0.6 0.8 1 0.2 0.4 Xi 0.6 0.8 1 




0.2 0.4 XiQ.6 0.8 1 0.2 0.4 2^^0.6 0.8 1 



Figure 3: Spatially non-uniform stationary solutions to (12.1 p which satisfy the conditions 
\xl — xj\ = 1, and the boundary conditions, (a) A solution that corresponds to the 
movement along half the cycle in Fig. (b) full cycle; (c) two full cycles; (d) four full 
cycles. Changing di and hence the velocity of the movement along the phase curves we 
can always obtain solutions with arbitrary number of full cycles 

then there exists a spatially non- uniform stationary solution to (12.11) . On the other hand, 
if d > ^ we obtain from Theorem 13. II that Ri are stable. Therefore we can consider d as a 
bifurcation parameter. As this parameter decreases spatially uniform stationary solutions 
become unstable, and spatially non-uniform solutions appear in the system according to 
the standard Turing bifurcation scenario. 

Now we consider the case of the spatially explicit hypercycle (12.21) . 

Theorem 4.2. Suppose that (14.11) holds. If the parameters of problem (12.21) can be rep- 
resented by one-parameter perturbation 

di = do -\- eli, = Oo + enii, rrii, li are constant, e > 0, 

where e is a small parameter, then there exist spatially non-uniform stationary solutions 
to system (12.21) . 

Proof. System (12. 2p can be rewritten in the following form: 

diu'- + UiioiUi - fi) = Uiai{ui - Ui_i) + (/i - /a)^,, i = l,...,n. (4.13) 

If e = then we have that /i = /a and 

dou- + UiiaoUi - /i) = 0, (4.14) 

which is a particular case of the autocatalytic system (12.70 . According to Theorem 14.11 
system (I4.14p possesses spatially non-uniform stationary solutions. Using the presentation 
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(I4.4p it can be shown that the right hand side of fl4.13p is of the order of e, i.e., can be 
rewritten in the form 

diu'- + Ui{aiUi - fi) = e'^iix), i = l,...,n. (4.15) 

where are bounded functions. This imphes that system fl4.15p is a perturbation of 

the Hamiltonian system (14.141) . According to the general theory [15] stable and unstable 
manifolds of the perturbed orbits will be close to the corresponding manifolds of the 
unperturbed system. Therefore for d < fi for each non-uniform stationary solution of 
(12. ip there exists spatially non-uniform stationary solution to (12. 2p . □ 



Remark 4.3. If we assume that the inverse to inequality (13.101) holds, then it can be 
rewritten in the form 

Indeed, we can rewrite inverse to (I3.10p in the form 



j=i 



1 

pn 



> 



1 

T 



Using the properties of arithmetic and geometric means we obtain 



n 

From the previous it follows that 



< n 



1=1 



ai p 



n 



i=l 



n 



Wd. 

i=l 



1 

pn 



n 

i=l 



n 



n 

i=l 



a,; 



> 



1 

pp 

2" 

TCP 



Once again using the inequality between arithmetic and geometric means we obtain 



n 



n 

i=l 



< 



i=l 



which proves the desired result. 

In words, we showed in this remark that if the inverse to (I3.10p holds, then the inverse 
to (14. ip is true, which means that if the spatially homogeneous solution to hypercycle 
system is stable there are no spatially non-homogeneous solutions. 
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Example 4.1. It is possible to obtain an explicit solution to (12.81) in the special case 
when n = 4, di = d^, ^2 = d^, ai = a, i = 1, 2, 3, 4. First, we rewrite fl2.11l) in the form 

/2 = / {Au,u)dx, U = (ui,U2,U3,U4y. 

Jo 

The expression {u,v) denotes the standard scalar product in M^, ' is the transformation. 
Matrix A is circular and has eigenvalues Ai = a, A3 = —a, A2 = A4 = 0. Consider the 
orthogonal transformation that reduces A to its canonical form: 



/ 1/2 


V2/2 


-1/2 


\ 


1/2 





1/2 


-V2/2 


1/2 


V2/2 


-1/2 





V 1/2 





1/2 


-V2/2 J 



Summing all equations in the hypercyclic system we have 

{uxx,D) = f2{u, 1) - {Au,u), 

where D is the diffusion vector, D = (^1,^2,(^1,(^2), 1 = (1,1,1,1)- Let u{x) = Tv{x), 
V = {vi,V2,V3,V4). It follows that 

{vx.,T'D) = f2{v,T'l) - {T'ATV,v). 

Since {T'ATV, v) = k{vf — fg), the last equation takes the form 

{di + d2){vi)xx = 2f2Vi - k{vi - v^). 

Suppose that M1+M3 = M2+M4. Then we have that f 3 = and the function w{x) = U1+U2, 
satisfies the differential equation 

{di + d2)wxx = 2/2W - aw'^, 

whose explicit solution can be found using (14. 4p . 

Remark 4.4. As in the case of system (12.11) parameter d can be considered as a bifur- 
cation parameter for (12.21) . 



5 Asymptotic behavior of the spatially explicit auto- 
catalytic and hypercyclic systems 

Consider the local system of autocatalytic reaction (II. ip in the form 

dwi 



dt 



^,(a,^f-/^(t)), /^(t) = ^a,^r\ t>s, 



i=l 



^i, z = l,...,n, ^Wi = l. 



i=l 

For the following we need 



(5.1) 
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Definition 5.1. We shall say that the initial conditions for system f l2.ip and system fl5.ll) 
are concerted if 



ii = Vi= I ^i{x)dx. 



(5.2) 



Let us assume that the initial conditions for systems (12. ip and (15. ip are concerted. On 
integrating system (12. ip with respect to x and using the equahty /S.v{x)dx = Jp d^vds = 
we obtain 

dvj 



dt 



ai I v^^^{x,t)dx-Vi{t)f[°\t), t>s, Vi{s) = ^i=^i, i = 1, 



cloc 



,n, 



where 



Viit) = / fi(x, t)dx. 
'a 



Since Ifil = 1 we have 



vf^'^{x,t)dx > I / Vi{x,t)dx 
n \Jn 



and, consequently, 

dv> 



>vf+\t)-v,it)f[°'it), t>s, v,{s) = ^, = ^,, t = l 



dt 



(5.3) 



Lemma 5.1. Let the initial conditions for systems (12.10 and (15. ip be concerted. Then 

< /{-(t) < f\{t), (5.4) 

_i 

where (3 = Ym=i " • 

Proof. First, we prove the left inequality in (15.40 . Using (12. 4p and Holder's inequality we 



have 



p+i 



< 



.1=1 



■ n ^ 

5Z 7^ 



i=l 



J2a.wr'=f3'fnt). 



i=l 



To prove the right inequality in (15.40 we assume that there exists s > that fi{s) < 
f["'^{s). Since the functions fi{t) and f[°'^ are continuous, there exists neighborhood 
Us = {t: < t - s < 6} from which /i(t) < /{°^(t) follows. Then from it follows 

that 

^>vr\t)-v,{t)f\{t), teUs, z = l,...,n. (5.5) 

Due to the fact that the initial conditions of (12.10 and (15. ip are concerted, then from the 
comparison theorem [28] we obtain 

Vi{t)>Wi{t), teUs, i = l,...,n, (5.6) 

where Wi{t) are the solutions to (15.10 . From the other hand we should have ^"^^^ Vi(t) = 
Sr=i ^j(^) = 1; we obtain a contradiction. 

□ 
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Theorem 5.1. Letp> 1. Then for almost all initial conditions (fi{x), Yll^=i lufii^)^^ ~ 
1 there exists an index j , I < j < n {which depends on (fi{x)) such that Vi{x,t) for 
all i ^ j in the space Lp^i, and J^Vj{x,t) 1 when t ^ oo. 

Proof. We have p > 1, and hence -ffp+i ^ The eigenfunctions ips, s = 0, 1, 2, . . . of 
the problem (13.21) form a complete system in H2. Let us represent 

00 

ipi{x) = c° + Zi{x), Zi{x) = y^^CiiJsjx). 

s=l 

Let Wi(t) be the solutions to (15. ip and let the initial conditions for systems (12. ip and (15.11) 
be concerted. We will look for a solution to (12.11) in the form 



(5.7) 



Vi{x, t) = Wi{t) + Zi{x, t), Zi{x, t) = ^ cl{t)ips{x), 

w,{0)=c^, cr(0) = cr, m = l,2,... 
Inserting (15.70 into (12. ID we obtain 

dwi{t) dzi{x,t) f u\/- /■4.\ I / I 7^^^j(^)^) 

-^jf- + — — = «»^f {x,t) - fi{t){wi{t) + Zi{x,t)) + di — — — . 

Integrating the last equation with respect to x and noting f^tps{x)dx = give 

= aj v^\x,t)dx - f,{t)w,{t). 
Using the fact that Wiif) are the solutions to (15. ip we obtain 

vt\xMx = iflit) - fT%t))w,{t) + a,wf^\t). {5.i 



It is known that solutions to (15. ip have a property of multistability. It means that all the 
vertexes of the simplex are stable, and the choice of initial conditions determines to which 
vertex the system evolves. In other words, for almost all initial conditions C,i the system 
(15. ip ends up in Rj, for which all the coordinates excluding j are zero {wi{t) when 
t — ^ 00 for all i 7^ j, and Wj(t) 1). Hence, from (15. 8p the theorem follows. □ 



Remark 5.1. Theorem (15.10 answers a natural question which spatially non-uniform 
stationary solution of (12. ip survives in the evolutionary process. To answer it we need to 
consider two systems (12. ip and (15. ip with concerted initial conditions. As was mentioned 
system (15. ip possesses the property of multistability; each vertex of the simplex has its 
own basin of attraction. If we denote these basins as Di, . . . ,Dn, then the number of 
the basin, to which the initial conditions of (15.10 belong, determines which spatially non- 
uniform solution will dominate the evolution. Note that for the dominant solution 

v{x,t)dx ^ 1 for t ^ 00. (5.9) 
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Another point here is that the exphcit space structure in the system with global 
regulation fl2.1l) does not provide the conditions for surviving more than one type of 
prebiotic replicators, Vi{x,t) — > in Lp+i for all i ^ j. 

Corollary 5.1. Almost all spatially non-uniform stationary solutions of the problem (12. ip 
are unstable. 

Proof. Consider (12. ip with the initial conditions 

Vi{x, s) = Ui{x) + Pi{x), i = 1, . . . ,n, 

where Ui{x) are spatially non-uniform stationary solutions to fl2.ip . and Pi{x) G Hs. From 
Theorem 15. II it follows that there exists a positive integer j (which depends on the initial 
conditions) such that Vi{x,t) ^ in space Lp+i for i j. Therefore only one set of 
stationary solutions can be stable. □ 



Remark 5.2. It is possible to obtain sufficient conditions for stability of the non- uniform 
stationary solution Uj{x) > 0, J^Uj{x)dx = 1. Unfortunately, applying this condition 
requires additional serious analysis. 

Indeed, we can look for a solution to (12. ip in the form 

Vj{x,t) = Uj{t) + Zj{x,t), Vi{x,t) = Zi{x,t), i ^ j. 

Putting these solutions into (12. ip and retaining only linear terms we obtain 

dtZj{x,t) = aj{p+ l)[u^{x)zj{x,t) — {Uj{x), Zj{x,t))] — ajZj{x,t){u^'^^{x), 1) + diAzj{x,t), 
dtZi{x,t) = -aiZi{x,t){u^i^^{x),l) + diAzi{x,t), i ^ j, 

with the initial conditions Zi{x,s) = Ui{x) + [3i{x). Here {u{x),v{x)) denotes the usual 
scalar product in L2(fi). This implies that all -2j(x, t) for i ^ j when t — > oo. On the 
other hand we have 

2di y^^i(^'^)^^ = l)[K(3^)'^j(3^'^)) - {uj{x),Zj{x,t)){u^j{x),Zj{x,t))] 

- aj{zj{x,t),Zj{x,t)){u^j^'^, 1) + di{Azj{x,t),Zj{x,t)). 
Substituting the following 

oo 

z,{x,t) = z';{t) + J2^jit)Mx) 

s=l 

into (IS.lOp and using the fact that 

oo 

{Azj{x,t),Zj{x,t)) = -Y^X^{z'.{t)Y 

s=l 
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we obtain that all the terms in (15.101) except for the terms in the square brackets are 
negative. The terms in the square brackets have the following form 



oo oo 



m=l s=l 



from which we obtain a sufficient condition for stability of the solution Uj{x) > in the 
form 



oo 



m=l 

The last formula should be checked only for small m because Xm ^ oo. 

The result of Lemma 15.11 can be extended to the case of hypercycle reaction. 

Lemma 5.2. Let the initial conditions of system ( \2.2\\ and system 



dt 



w,{aiwl, - ftit)), ftit) = J2 t > 0, 



1=1 



u!i{0)=^i, i = l,...,n, y^Wi = 1 



(5.11) 



i=l 

be concerted. Then 



ftit) < /2(t). (5.12) 



Proof. We have 



and 



AVi f dVi 1 I 1^ 1 ly 7 

dx = / -;- ds + I \\/xi^Vi\ dx > 0, 

^ Vi Jr on Vi 



Therefore 



vf_-j^{x,t)dx > Vi^i{x,t)dx^ = '^t 
f d 

/ — In Vidx> aivf_-^{t) - f2{t). 

Since the initial conditions of (12. 2p and (15. lip are concerted, then, as in the case of 
Theorem (15. ip we can represent fi(x, t) as the sum fj(x, t) = Wi(t) + Zi{x, t), where Zi{x, t) 
are given by (15. 7p . and note that WiifS) = Vi{0) for any i. 
From the last inequality it follows that 

d d _ 

— Invi dx = -^^Wiit) > aiV^_^ - f2(t). 

Since Vi{0) = Wi{0) then, using (15.111) we obtain (I5.12p . □ 
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Using the last lemma we can extend the results of permanence of hypercycle system 
with p = 1 to the spatially explicit case [18j. We remind that permanence means that 
solutions to system (15.111) with the initial conditions Wj(0) = .^j > do not vanish, i.e., 

1 > Wiit) > 5 > 0, t > 0. 

Corollary 5.2. Let p = 1 and let the initial conditions of systems (12.21) and (15.111) be 

concerted, and 

^i = Ci = MO)>0, J]ei = i- 
Then the solutions to system ( 12.21) do not vanish in L2 space. 
Proof. Let a solution Vi{x,t) to (12.21) vanish for some i, i.e., 

\\Vi{x,t)\\L2 0, t 00. 



Using the reasoning along the lines of Theorem 15. H we obtain 

Vi{x,t)vi-i{x,t) dx = (/2(t) - f2''{t))wi{t) + aiWi{t)wi^i{t). 
Using Lemma [5.21 we hence have 

Vi{x,t)vi-.i{x,t) dx > Wi{t)wi^i{t) . 
The last and the Cauchy inequalities yield 

\\Vi{x,t)\\L2\\Vi-iix,t)\\L2 > Wi{t)Wi_i. 

From the fact ||fi(a;, t) — > it follows that either Wi or Wi-i tend to zero, which 
contradicts to the permanence of the hypercycle system (15.111) . This completes the proof. 

□ 

Similar to Remark 5.2 we can obtain sufficient conditions for stability of the spatially 
nonhomogeneous stationary solutions for the hypercycle system (12.21) . However, the utility 
of such conditions is questionable because we hardly can expect that we will be able to 
check these conditions analytically. 

It is possible to study the stability of spatially nonhomogeneous solutions in somewhat 
weaker sense. 

Definition 5.2. We shall say that spatially non-uniform stationary solution u{x) = 
{ui{x) , . . . , Un{x)) to system ( 12. ip or ( 12.21) is stable in the sense of the mean integral 
value if for any e > there exists 6 > such that for the initial conditions 

\ipi - Ui\ < 6, 

it follows that 

\vi - Ui\ < e, 

for any i and t > 0, where, as before, Vi{x,t) are the solutions of (12.11) or ( 12.21) . 

Vi = Vi{x,t)dx, Lpi = I ipi{x)dx, Ui = Ui{x)dx. 



20 



It is clear that the stabihty in the mean integral sense is weaker than the stability in the 
usual sense (Lyapunov stability). For example, consider functions g{x,t) G H^, x G [0, 1] 

oo 

g{x, t) = co(t) + Cfe(t) cos knx. 

s=l 

Let us suppose that Co(t) — ^ when t — > oo. Then g(t) whereas t)||^i = 

Yl"^! + fc^TT^) does not necessarily tend to zero. 

Corollary 5.3. Let us suppose that the following inequalities hold for any i = 1, . . . ,n: 

*<f, (5.13) 

at Ai 

Then all spatially non-uniform stationary solutions to fl2.ll) of the form 

f/^(x) = (0,...,0,M,(x),0,...,0) 
are stable in the mean integral sense. 

Proof. From Theorem 13. II it follows that Rj are unstable when fl5.13p holds. Consider the 
solution W{x) for which uj = 1. From the other hand from Theorem (15. ip follows (15.90 . 
which completes the proof. □ 

Now we switch to the hypercycle system (12. 2p with explicit spatial structure and global 
regulation. After integrating (12. 2p with respect to spatial variable, we obtain 

^ = ai{vl„ V,) - f2{tMt), < t, v,{0) = (5.14) 

where the meaning of the function g as before, in the mean integral sense. 
Let us introduce new functions 

Vi{x, t) = i = l,...,n, 0= / y ^{aj)pVj{x,t)dx. (5.15) 

For the new variables 



= / y ^{aj)pVj{x, t)da 



n „ 

y / Wj{x,t)dx 

,=1 



(5.16) 



Note that in the new variables the stationary point P2 has the coordinates (1/n, 1/n, . . . , l/ra) . 

Lemma 5.3. In the new variables (I5.15P the dynamical system (I5.14p has the following 
form: 

^ = ent)(K-„i.,) -/2(t)^.(t)), 

" f (5.17) 
h{t) = 2. / Wj{x,t)w^_^{x,t)dx. 
3=1 
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Proof. Using (15.151) . (12.21) and the boundary conditions we obtain 



dQ 
'dt 



n „ 

n „ n 



Wjdx 



Using (I5.16P we obtain 



de 

'dt 



n „ 

6^+^ V / Wjw'l_^dx - /2(t)e. (5.18) 



Equality (I5.15P yields 



dvi WiQ + WiQ 

From the other hand (15.151) implies 

dVj ^ QP+\Wi,wU)-QWrUt) 



(5.19) 



(5.20) 



Putting together (I5.18p . (l5.19l) and (I5.20p completes the proof. □ 

Consider the spatially uniform stationary solution P2 to (12. 2p . It is also an equilibrium 
of fICTD . 

Theorem 5.2. Let u{x) = {ui{x), . . . ,m„(x)) be a spatially non-uniform stationary solu- 
tion to (12. 2p such that u = P2, where P2 is the homogeneous stationary solution of (12.20 . 
Then u{x) is stable in the sense of the mean integral value. 

Proof. Consider system (I5.14p . Due to Lemma [5. 31 this system is topologically equivalent 
to system (I5.17p . which has the steady state Pq = (1/n, . . . , l/ra). Let us introduce the 
following Laypunov function 

V{wi, . . . , Wn) = — ln{wiW2 . . . Wn) — nluTi. 

It is easy to see that V{Pq) = and V{wi, . . . , Wn) > in a neighborhood Zs of Pq, where 



Wi,i = 1, . . . ,n 



n " 1 1 
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Using fl5.14p yields 



V 



i=l 



i=l 



{Wi, wf_^) 



Wi 



n 

-e^5^(«;.,t.f_,)(--n). 

. -I Wi 

1=1 



Denote the following 



/i = min |inf(wi,«;f_i)| . 
The functions Wi{x,t) are nonnegative for all i, therefore we obtain 



We also have 7^7-1 ^ > 



Since X]r=i = 1; "^^i ^ 0, the function HILi ^« 



reaches its maximum at the point Pq = {1/n, . . . ,1/n), and this implies that 

n 



Wi 



which means that < 0. Invoking the arguments of the topological equivalence of fl5.14l) 
and fl5.17p completes the proof. □ 



6 Conclusion 

In this paper we studied the existence and stability of stationary solutions to autocatalytic 
and hypercyclic systems (12.1 p and (12. 2p with nonlinear growth rates and explicit spatial 
structure. It is well known that the mean field models (e.g., models described by ODE 
systems) are often show different behavior from the models where the spatial structure 
is taken into consideration (more on this [8j). In particular, it is widely acknowledged 
that the evolution and survival of altruistic traits can be mediated by spatial heterogene- 
ity. Macromolecules that catalyze the production of other macromolecules are obviously 
altruists, and in this note we tried to answer the question whether the particular form 
of spatial regulation (namely, global regulation [H [32]) can promote the coexistence of 
different types of macromolecules in the prebiotic world (within a hydrothermally formed 
system of continuous iron-sulfide compartments [21] )• The analysis presented in [3ii32j is 
significantly extended to the cases of nonlinear growth rates, arbitrary fitness and diffusion 
coefficients. 
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0.6 

X 0.4 



Figure 4: The competitive exclusion for autocatalytic growth. Numerical solutions to 
autocatalytic system (12. ip . n = 3, rfi = 0.02, d2 = 0.05, = 0.08, p = 1, ai = a2 = 
= 1. The initial conditions are Ui{x,0) = 0.35 + 0.3cos7ra;, ^2(0;, 0) = 0.35, Us{x,0) = 
0.3 — 0.25 cos vrx. Note that the orientation of the axis is different for (a) and (6), (c). 
Only one type, ui, survives. The asymptotic state is a spatially non-uniform stationary 
solution. The details of the numerical computations are given in [1] 



The major conclusion is as follows: the mathematical models with spatial structure 
and global regulation show in general very similar qualitative features to those of local 
models. Two basic properties, namely the competitive exclusion for autocatalytic systems 
and the permanence for the hypercyclic systems, are shown to hold for spatially explicit 
systems. Numerical calculations illustrate these conclusions in Figs. H] and [5] (the details 
on the numerical scheme used in the calculations are given in [1]). 

More precisely, for sufficiently large diffusion coefficients the spatially uniform station- 
ary solutions to (12. ip and (12. 2p have the same character as in the local models (II. ip and 
(II. 2p . For such diffusion coefficients the asymptotic behavior of the local and distributed 
models coincides. If, on the other hand, the inequality (14. ip holds and the nonlinear 
growth rates satisfy the condition < p < 2 then new, spatially non-uniform solutions 
appear; for small diffusion coefficients these spatially heterogeneous solutions can corre- 
spond to the multiple cycles on the phase plane of the corresponding Hamiltonian system 
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1 



Figure 5: The permanence for the hypercycle system. Numerical solutions to hypercyclic 
system (Q- n = 3, = 0.001, rfa = 0.002, = 0.003, p = 1, ai = as = 03 = 1. The 
initial conditions are Mi(x,0) = 0.35 + 0.15 cos vrx, ^2(0;, 0) = 0.357, M3(x,0) = 0.338 — 
0.3 cosvra;. The asymptotic state is spatially non-uniform stationary solutions. The details 
of the numerical computations are given in [1] 



(Fig. [3]). In the case of autocatalytic system these solution can be stable only if all but 
one asymptotic state are zero. In the case of the hypercyclic system we prove that these 
spatially heterogeneous solutions can be stable in the sense of the mean integral value. 
The examples of the asymptotic states for a hypercyclic systems found numerically are 
shown in Fig. El These non-uniform stationary solutions can be considered as the means 
of the hypercycle system to withstand the parasite invasion [22] (the analysis of models 
with parasites and with p > 2 is the subject of the ongoing work). 

A Appendix 



Proof of Lemma \4-l\ Consider the function g{t) = — p > 0. This function has 



two roots Ti = and r2 = 1, and attains its maximum at t* 



2+p 



, which is g{t*] 
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0.2 0.4 X 0.6 0.8 1 0.2 0.4 X 0.6 0.8 1 




Figure 6: Asymptotic spatially heterogeneous states of the hypercycle system f l2.8p found 
numerically, n = 3, di = 0.001, d2 = 0.002, = 0.003, p = 1, ai = 02 = as = 1. Note 
that case (/) corresponds to the simulation shown in Fig. [5l 



2+p 



. Function Pi{t) can be obtained from g{t) by shifting the latter. Therefore, 
when (14. 5 p holds, Pi{t) has two positive roots that are situated in the interval (14.60 . □ 

Proof of Lemma \4-^ To simplify notations we drop indexes where it is possible. We need 
to prove that for 

P(r, c) = c + T^- tP+^ 



and 

h{c) - 
where 

Ti G (0,ro), T2 e (to, 1), To 
we have that 



dr 



h{c) 



r, c 



rdr 



(A.l) 



2+p 



, P{n,c) = 0, P{t2,c) = 0, P;(ro,c) = 0, 



/(c) = -{hic))l-'hic) > 

To 



P 



(A.2) 
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for < p < 2. 

For p — 2 direct calculations show that /2(c) = |, /(c) = hence we assume that 
< p < 2. Using Holder's inequahty yields 



where 



h(c) 



T2 



dT 



c) 



dr 



Next we will the following change of the variables: 



E+2 

r 2 



P+2 



2±2 



p + 2 



g(i,c) = c+goW, Qo{t)^t''-t\ i<?<2, 



from which c) = Q{t2, c) = 0, and hence c = — (5o(^i) = "^0(^2)- 
Out integral takes the form 



h(c) 



t2 



dt, 



tie{0,to), t2e{to,i), to 



q\2-- 



, Qo(to) = 0. 



Function Q(t, c) does not exceed its Hcrmitc interpolation polynomial H^, which is 
build using the values at H'^{ti) = H^{t2) = H'^{to) = 0, H-i{tQ) — Q{to,c). This follows 
from non-negativity of the reminder term of interpolation 



Q{t,c)-Hs{t)^ 



and the fact that QI^\0 > when ti < ^ < t2. Therefore, we have 
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{t-to)\t-ti){t-t2), 



dt, 



where 



Hs{t)^Q{to,c)ll 



it-to){2to-ti-t2)\ {t-tl){t-t2) 



{to - tl){to - t2) J {to-ti){to-t2)' 

Making the change of the variable in the integral 

ti + t2 t2 — ti . 

t — — 1 — sm ip, 



we obtain 



t2 



dt 



{to-ti){t2-to) 



Q{to,c) 



h{c), 
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where 



J-n/2 Jl-^ 



+t2)/2-f +(t2 -ti )/2 sin y)(2to -ti -ta) 



{tQ-tl){tQ-t2) 

Since the graph of any convex function lays above any tangent hne, then we have 

1 ^ / CiX 



VClX + C2 V 2C2 

for any x. Using the last inequality we can estimate I 4 as 

I4 > / ^ 1 d(p = —— = > TT. 

J-KI2 \/C2 V 2C2 y /. (2to-ti-t2)^ 

V 2(to-ti)(t2-to) 

Using the last estimate and returning to I3 we obtain that 



r ^ \ ^ 1^ {to - ti){t2 - to) qn ^ . , ^ . , x 



where 



Qo(io) - Qo{t) ■ 

With the help of the Taylor formula the denominator in g{t) can be presented in the 
following form: 

Qo{to) - Qo{t) = -Q'^ito)it - to)V2 - Qoito)it - tof/Q - Q?^(C)(t - to)V24, 
where ( belongs to the interval {t, tg). If we denote C3 — Q'o{to)/6, we obtain 

^^^'^^^^'^ ^ (2-g-C3(ti-to))(2-g-C3(t2-to)) 
Denominator of this fraction 

(2 - qf + C3(2 - q){2to - h - ^2) + cl(ii - to){t2 - h) 

has its fist term positive and its second and third terms negative. Indeed, we have 
Qo(^i) = Qo(^2), and, using the Taylor formula around t — to for both parts of this 
equality, we obtain 

(g - 2)(ti - to)' + Qo'(ei)(ti - to)V6 = {q- 2)(t2 - to)' + Qo"(6)(t2 - tof/6, 
where ^1 e {ti,to) and ^2 G (^0)^2)- Then 

(g - 2)((ii - io)' - (^2 - toY) = Q'^{m2 - tof/G - g'o"(6)(ti - to)V6 < 0, 
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since Q'^if) < for any t. Which imphes that (ti — t2)(^i + ^2 — '^io) from which follows 
that the second term is negative. Using this fact we obtain 



/ X / X 1 T Q 



^0 ^0 V P / VP + 2^ 

which completes the proof. □ 
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